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INTRODUCTION 


This report describes the current status of work at Wallops Flight Center 
In the past year during which over-ocean data from the waveform samplers In the 
SEASAT-1 radar altimeter have been carefully analyzed In attempts to extract the 
skewness of the ocean surface elevation distribution. As described later, our 
work also provides Information on possible tracking-point shifts due to changes 
In the surface distribution; a track-point shift Is a type of altitude measure- 
ment error. The basic description of the radar altimeter Is provided by Ref. 1, 
and Ref. 2 summarizes Its performance evaluation. Preliminary results from our 
work had been contributed to the GOASEX-II Workshop In June 1979 (Ref. 3) and to 
a SEASAT Height/Orbit Accuracy Workshop In November 1979. Since then, we have 
spent considerable time trying to get around the difficulties of the combination 
of poor knowledge of the off-nadir attitude angle and the Individual waveform 
samplers gain uncertainties, problems which will be described more fully later 
In this report. 

There were several different motivations for this work. First, the problems 
we have encountered in the SEASAT-1 data analysis should be of immediate concern 
to the radar altimeters to be carried by the NOSS oceanographic satellite, and 
the data processing procedures we develop will carry over to the NOSS system. 
Second, one of the largest uncertainties in the SEASAT Instrument height accuracy 
is the uncertainty in the possible track-point shift in response to changes In 
the ocean surface elevation distribution. Third, the surface skewness Is a data 
product of considerable interest to oceanographers since, in combination with 
the significant waveheight (SWH), it allows estimating the ocean surface domi- 
nant wavelength based on the work of Huang and Long CRef. 4). (See Ref. 3 or 
Ref. 5 for examples.) We are in the process of working on data from some of the 


storms listed in the recent SEASAT Storms Workshop report (Ref. 6) and will 
report our results in a later publication. 

The work reported here has been performed in collaboration with D. W. 
Hancock, III and with continuing help from E. J. Walsh. Additional assistance 
has been provided to us by R. G. Forsythe, J. Newell, and G. Quintana, and we 
have benefited from useful discussions with N. E. Huang, W. F. Townsend, C. L. 
Purdy, J. T. McGoogan, and J. L. MacArthur. The background for some of the 
waveform fitting lies in work performed during the research on SWH extraction in 
the GEOS-3 project under H. R. Stanley. 

In the rest of this report we discuss the basic approach used, the waveform 
fitting procedure, and the model waveform fitted based on a near-Gaussian point- 
target response. The SEASAT-1 actual point-target response and its consequences 
are then discussed, followed by waveform sampler gains and attitude angle 
effects. Finally a brief summary section concludes the report. 

BASIC APPROACH IN WALLOPS WAVEFORM ANALYSIS 

The SEASAT-1 radar altimeter has a set of 63 waveform samplers which pro- 
vide, once every tenth of a second, an average of the last 100 individual radar 
return waveforms. Sixty of these samplers are uniformly spaced, providing 
samples at every 3.125 nanoseconds along the radar return waveform, and the 
remaining three samplers are spaced 3.125 ns from each other and lie at the 
center of the set of 60 to increase the sampling density at the most important 
part of the return waveform (the "ramp" region). Table 1 lists the time loca- 
tion of these 63 samplers and establishes the numbering system for referring to 
them; we will refer to samplers 1 through 60 plus numbers Z9it 30i, and 3H. In 
addition to the pre-launch data, two types of on-orbit calibration data are 
available for the 63 waveform samplers. In Calibration Mode I, the transmitted 


TABLE 1 . TIME LOCATION AND INDEXING FOR THE 63 SEASAT WAVEFORM SAMPLERS 


Previous #, This 


Index 

Time, ns 

SEASAT # 

Report 

1 

-92el875 

-30 

1 

2 

-89.0625 

-29 

2 

3 

-85.9375 

■28 

3 

4 

-82.8125 

-27 

4 

5 

-79.6875 

-26 

5 

6 

-76.5625 

-25 

6 

7 

-73.4375 

-24 

7 

8 

-70.3125 

-23 

8 

9 

-67.1875 

-22 

9 

10 

-64.0625 

-21 

10 

11 

-60.9375 

-20 

11 

12 

-57.8125 

-19 

12 

13 

-54.6875 

-18 

13 

14 

-51.5625 

-17 

14 

15 

-48.4375 

-16 

15 

16 

-45.3125 

-15 

16 

17 

-42.1875 

-14 

17 

18 

-39.0625 

-13 

18 

19 

-35.9375 

-12 

19 

20 

-32.8125 

-11 

20 

21 

-29.6875 

-10 

21 

22 

-26.5625 

- 9 

22 

23 

-23.4375 

- 8 

23 

24 

-20.3125 

- 7 

24 

25 

-17.1875 

- 6 

25 

26 

-14.0625 

- 5 

26 

27 

-10.9375 

- 4 

27 

28 

- 7.8125 

- 3 

28 

29 

- 4.6875 

- 2 

29 

30 

- 3.1250 

- H 

29i 

31 

- 1.5625 

- 1 

30 

32 

0.0000 

0 

— 

o 

33 

1.5625 

+ 1 

31 

34 

3.1250 

+ H 

3H 


Index 

Time, ns 

Previous 
SEASAT # 

#. This 
Report 

35 

4.6875 

+ 2 

32 

36 

7.8125 

+ 3 

33 

37 

10.9375 

+ 4 

34 

38 

14.0625 

+ 5 

35 

39 

17.1875 

+ 6 

36 

40 

20.3125 

+ 7 

37 

41 

23.4375 

+ 8 

38 

42 

26.5625 

+ 9 

39 

43 

29.6875 

+10 

40 

44 

32.8125 

+11 

41 

45 

35.9375 

+12 

42 

46 

39.0625 

+13 

43 

47 

42.1875 

+14 

44 

48 

45.3125 

+15 

45 

49 

48.4375 

+16 

46 

50 

51.5625 

+17 

47 

51 

54.6875 

+18 

48 

52 

57.8125 

+19 

49 

53 

60.9375 

+20 

50 

54 

64.0625 

+21 

51 

55 

67.1875 

+22 

52 

56 

70.3125 

+23 

53 

57 

73.4375 

+24 

54 

58 

76.5625 

+25 

55 

59 

79.6875 

+26 

56 

60 

82.8125 

+27 

57 

61 

85.9375 

+28 

58 

62 

89.0625 

+29 

59 

63 

92.1875 

+30 

60 
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radar pulse Is fed directly to the sampler set through a variable attenuator; 
this provides the best available knowledge of the system point-target response 
function. In Calibration Node II, a uniform signal Is presented to the sampler 
set. Each sampler has a multiplicative gain differing slightly from Its neigh- 
bors and Calibration Mode II Is Intended to provide the necessary gain correc- 
tion Information. Our work Is based on the notion of using data from all 63 
samplers, averaged over some suitable smoothing Interval, and "best-fitting" 
these 63 experimental values to a model waveform by varying the model waveform's 
parameters until a minimum Is reached In the sum of squares of the llnefit 
residuals. Our usual smoothing Interval Is about 10 seconds, so that at SEASAT's 
rate of 1000 pulses per second, a waveform is fitted to a 10000 pulse estimate 
of the mean return. 

The success or failure of this wavefornsfittlng will depend primarily upon 
the adequacy of the model waveform at representing the actual sampled waveform 
data being fitted, assuming that: 1) the parameter-varying, the fitting itself, 
is performed adequately, and i1) that the correct gain calibrations are known 
and have been applied. Our work uses a model waveform described by seven 
parameters: i) amplitude; ii) time origin ("track-point"); iii) composite, or 
total, risetims; iv) baseline; v) composite, or total, skewness; vi) off-nadir 
attitude angle; and vii) composite, or total, kurtosis. In all work to date we 
have used a fixed zero kurtosis. 

The six parameters fitted are necessary descriptors of the system-observed 
mean return waveform. The baseline for example is obviously present in the data 
and so must be fitted for, but is of little interest in itself. Listed below 
are several of the fitted parameters with comments in their possible uses or 
implications. 
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Amplitude - Not of much interest for now, but should be included in any 
attempts to refine the ocean surface reflectivity (c®) measurements. 
Attitude Angle - Of interest because of its use in estimating a®. One can 
derive a theoretical a® correction as a function of attitude. It was 
anticipation of the a®-related need for good attitude estimates which led 
us originally to try to use all 63 waveform samplers in the fittingj more 
samples used from later in the vmveform plateau should lead to more stable 
estimates of attitude. 

Composite Risetime - This is used in the SWH estimation. The SEASAT-1 on- 
board processor produces estimates of SWH which are probably better than 
any user's current ability to use SWH; nonetheless, the six-parameter wave- 
form fitting will lead to a slightly different, somewhat more accurate SWH 
value than the on-board processor. 

Composite Skewness - This contains the ocean surface skewness, an oceano- 
graphic parameter of considerable potential interest. We do not yet know 
the limiting stability of the ocean surface skewness estimated from the 
SEASAT-1 experimental data but will be determining the skewness stability 
in our continuing waveform analysis research. 

Time Origin - This parameter indicates the fitted waveform position rela- 
tive to the set of waveform samplers, and hence can be used to provide a 
correction to the altimeter's altitude output. The waveform fitting pro- 
cedure can be viewed as providing a "software tracker" which replaces, 
after the fact, the results from the SEASAT-1 hardware tracker. Possible 
shifts in the hardware tracking point due to changes in the ocean surface's 
distribution will show up as changes in the time-origin parameter. 
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WAVEFORH FITTING ROUTINE 


The model waveform to be fitted Is supplied as a subroutine to the overall 
fitting procedure. The entire fitting routine Is quite general and permits: 1) 
varying the number of experimental data points being fitted; 11) varying the 
number and order of the waveform parameters being fitted; 111) applying con- 
straints to one or more of the fitted parameters; and 1v) uniform or non-uniform 
weighting of the input data. In the work to date* the experimental data have 
been uniformly weighted. We have found that the solutions seem more stable If 
the skewness and the pointing angle are slightly constrained when they are among 
the fitted parameters. The overall fitting routine Is a straightforward generali- 
zation and extension of the approach sketched In Ref. 7 for a four-parameter 
situation. A description of the general fitting program will be published as a 
separate report, but a Fortran source listing and further comments may be 
obtained In the meantime by contacting either G. S. Hayne or D. W. Hancock, III 
at NASA Wallops Flight Center. 


MODEL WAVEFORM FinED 

Based on earlier work of a. number of researchers as summarized In Ref. 8, 
It Is assumed that the correct mean return waveform from a typical radar altim- 
eter can be obtained from the convolution of the following three terns: 1) the 
average Impulse response function of the quasi-calm ocean surface (the "flat 
sea"), 11) the ocean surface elevation distribution function for the "effective 
radar scattering elements," and ill) the radar system point-target response 
function. The first of these terms, the quasi-calm ocean impulse response 
function, is provided by Ref. 8. It looks somewhat like a unit step function 
with a subsequent "plateau" behavior which is a function of the antenna beam- 
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width and off-nadir angle. Figure 1 shows the flat-sea Impulse response func- 
tions for several different off-nadir angles for SEASAT-1 and for GEOS-3. The 
location of the "altitude gate," a wide waveform sampling gate, Is shown for 
later discussion. Note that the Figure 1 response functions have been renor- 
malized to Include (approximately) the effects of the SEASAT-1 and GEOS-3 AGC 
systems. In this report we are concerned with the shape of the Impulse response 
functions end the waveforms, not the amplitude as such. The surface reflec- 
tivity, a®. Is treated In the rest of this report as an arbitrary unknown multi- 
plier; Its value can be estimated through analysis of the AGC signal, but this 
report will not treat the a® estimation further. 

The second term In the convolution Is the radar-observed surface elevation 
probability density function (pdf) and In past work It has been assumed that the 
radar-observed elevation pdf Is the same as the true, geometrical surface eleva- 
tion pdf. I have discussed this briefly In Ref. 9. Ref. 10 Is a recent paper 
which treats the differences In distribution on a theoretical basis, but for the 
work In this report it Is sufficient to assume that there is some distribution 
varying with SWH and to put off these more serious questions of what distribu- 
tion until later. 

The third term in the convolution, the system point-target response, is the 
radar altimeter's transmitted pulse shape as sampled at the receiver by the 63 
waveform samplers in the SEASAT Instrument. We use the phrase "system point- 
target response" here instead of the more casually used "radar pulse shape'^' 
because what is meant is the transmitted radar pulse as modified by receiver 
effects . The reader might for ease of physical interpretation prefer to think 
of the point-target response as the same as the transmitted pulse, and this is 
nearly correct. 


Figure 1. Calm sea Impulse response functions vs time for GEOS-3 
and SEASAT-1 radar altimeters. (These functions are 
zero for all times less than zero.) 
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The convolution of the three terms Just described Is shown schematictslly In 
Figure 2 In which Pp 5 (t) Is the function alr*^a# shown In Figure 1. Notice that 
all three terms In the convolution are assumed to bu time functions. The 
surface elevation pdf Is a spatial distribution but can be transformed to a time 
function by using the altimeter's (two-way) ranging velocity of light* c/2 = 0.15 
meters per nanosecond. In Figure 2. both the surface elevation pdf* q^(t) and 
the radar system point-target response* s^(t)* are sketched as being nearly 
Gaussian. In this report we will refer to a function's being nearly Gaussian as 
meaning specifically that the function* q^(t) for example* can be adequately 
represented by the skewed Gaussian form 

q^lt) - exp [- (t/Oj)^] . 

where Is the second moment (the "width") and Is the skewness. The skew- 
ness Is dimensionless and has time units, nanoseconds In this work. If the 
system point-target response function s^(t) Is also nearly Gaussian with width 
and skewness X^, then the convolution s^(t)*qg(t) will also have this same 
nearly Gaussian form (neglecting higher moments) with a composite width and 
composite skewness X given by 

a = [a ^ + a 
c *■ s r 


and 


X 


c 
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A separate paper (Ref. 9» also Ref. 10) gives a several-tem series expan- 
sion which for GFASAT-1 adequately approximates the convolution of a nearly 
Gau^;s1an form with the Pp 3 (t) of Figure 1. The advantage of the results In 
Ref. 9 Is that the twc^ convolutions Indicated In Figure 2 can be replaced by 
expressions whose evaluation requires much less computer time than would numeri- 
cal convolutions. 

Figures 3 and 4 are reprinted from Ref. 10 and show the dependence of the 
mean return waveform on SWH and on attitude angle for a radar altimeter having a 
pure Gaussian point-target response function with the SEASAT-1 nominal design 
pulsewldth, and with the SEASAT-1 antenna beamwidth and satellite altitude. 
Figure 5, also reprinted from Ref. 10» shows the difference that a surface 
skewness will make In the radar altimeter's mean return waveform. A skewness 
magnitude of 0.5 Is about the outside limit that anyone might expect for a real 
ocean surface elevation distribution* so Figure b shows the outside limits for 
the skewness-caused waveform changes In the SEASAT case and shows something of 
the difficulty of any attempts to measure a surface skewness. 

All of the Wallops waveform work Is based on the wavefom; expansion given 
In Ref. 9. In the event that the radar system point target response function 
cannot be adequately represented by a nearly Gaussian function* we still use 
the expressions given In Ref. 9 to represent the convolution of the flat-sea 
impulse response with the surface elevation pdf (convolution of Pp 5 (t) and qg(t) 
In Figure 2)-, a numerical convolution of this result with the point target 
response (s^(t) In Figure 2) then gives the complete waveform. In other words* 
we use the Ref. 9 expressions to represent "everything else*" everything but the 
point target response function* and then convolve the "everything else" function 
with the point-target response for the complete waveform. With the Ref. 9 
waveform as the starting point for the Wallops work, the difficulties encoun- 


Figure 3. SEASAT-1 mean return waveform change with significant 
waveheight. 
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MEAN RETURN WAVEFORMS & DIFFERENCES 


Figure 5. SEASAT-1 mean return waveform change with skewness. 
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tered Included 1) the actual shape of the po1nt>target response, 11) Individual 
sampler gain biases, and 111) knowledge of the off-nadir attitude angle. 

These Items will be described In following sections of this report. The terms 
“off-nadir angle," “attitude angle," or “attitude" are treated as Interchange- 
ably synot\ymous with each other and with “altimeter pointing angle off-nadir" In 
the following discussions. 

ACTUAL SEASAT-1 POIlff-TARGET RESPONSE AND ITS CONSEQUENCES 

The SEASAT-1 nominal design system point-target response was 3.125 nano- 
seconds wide at 1/2-power points, and the on-board processor for SWH was de- 
signed assuming this to be a Gaussian shape with a second moment value 1.327 ns 
(the ratio of the half-power full -width to the second moment of a Gaussian 
function Is 2.355; 3.125 ns t 2.355 = 1.327 ns). The unweighted chirp trans- 
mitted radar pulse Is supposed to have a [(sin x)/x] form, and a first-order 
correction to the SWH to account for the non-Gausslan shape Is described In 
Ref. 1. This correction was simply to replace the 3.125 ns pulse by a Gaussian 
function 20% wider In the SWH calculation, and the SEASAT data processing at JPL 
Includes a look-up table which makes this correction. Figure 6 shows the actual 
system point-target response averaged from eight sets of Calibration Mode I, 

Step 9 data from the SEASAT-1 waveform samplers. The only half- Integer-numbered 
waveform samplers are numbers 29i, 301, and 311 (see Table 1 for notation), and 
In Figure 6 we have put zero values (except for the noise baseline) at every 
other half-integer location (where nulls should have occurred In the (sin x)/x 
pulse). This figure, assuming the Individual sampler gains are correct, repre- 
sents the best available knowledge of the actual system point- target response 
function. Because the point-target response In Figure 6 Is obviously not 
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Figure 6. SEASAT Radar Altimeter System Point-Target Response 
from CAL Step 9 after Gain Bias Corrections 
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Gaussian, not [(sin x)/x] , not even symmetric, a waveform fitting procedure 
based on a near-Gaussian point- target response will be fitting the data to the 
wrong model waveform with the consequence that every fitted parameter will be in 
error. These errors decrease as SUM increases (as the surface distribution's 
second moment becomes appreciably larger than the second moment of the point- 
target response). Figures 7 and 8 show the results of fitting the Gaussian- 
based model waveform of Ref. 9 to the mean return waveform samples which would 
be the result of a near-Gaussian surface distribution and the system point- 
target response shown in Figure 6. The "data," the individual points in Figures 
7 and 8, are generated by a numerical convolution of the point-target response 
of Figure 6 with an "everything els?" function from Ref. 9, and the solid lines 
in Figures 7 and 8 are the fitted waveforms using the incorrect model waveform 
(that is, the Ref. 9 results based on a near-Gaussian point-target response). 
Every parameter fitted has some degree of error, and this is particularly 
extreme in Figure 7 for zero input SWH; the incorrect solution gives an off- 
nadir angle of nearly 0.4" although the true input angle was 0.0". 

The true point- target response of Figure 6 can be built into the line- 
fitting procedure by incorporating within the model waveform subroutine a 
numerical convolution of the point-target response and the "everything else" 
functions; the penalty for this is considerably longer computation times in the 
waveform analysis programs. Our waveform analysis programs were written for 
research purposes with a high degree of flexibility and for ease of model func- 
tion modification and no particular effort has been directed at writing effi- 
cient source code for reducing computation times, but some idea of the time- 
penalty of the numerical convolution within the model waveform can be gotten 
from our program's running times on a medium-sized digital computer. A six- 
parameter fit of the near-Gaussian point-target response function model waveform 
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(represented solely by the expressions of Ref. 9) to 63 waveform samples ran In 
somewhat less than one second. In contrast, the "correct model" Incorporating 
numerical convolution with the point target response of Figure 6 required 30 
seconds or more to fit the same data. 

This discussion has been concerned with two different considerations: 1) 
adequacy of representing the actual expected mean wavefomu and 11) computer 
running time. We Indicated two different approaches which. In the rest of this 
report, will be referred to as 1) the "Incorrect but fast" model waveform, and 
11) the "correct but slow" model waveform. The "Incorrect but fast" waveform 
assumes the point-target response to be Gaussian and then uses the general 
expressions of Ref. 9 to give the waveform directly. The "correct but slow" 
waveform uses the Ref. 9 expressions as the "everything else" function which Is 
then numerically convolved with the actual point- target response of Figure 6 to 
give the final waveform. 

CORRECTION CURVES ALLOWING USE OF "INCORRECT BUT FAST" MODEL WAVEFORM 

This paragraph will describe one possible way to avoid paying the penalty 
of longer computation times and still take account of the effects of the real 
system point-target response function. It should be possible to use the "incor- 
rect but fast" waveform model (as defined in the last paragraph above) in the 
waveform fitting and then correct the resulting fitted parameters by using a 
separately-generated correction table. This possibility is suggested by Figures 
7 and 8 on which the correct (input) parameters and the incorrect (fitted) 
values are listed; for the particular cases in these figures we have the neces- 
sary correction data. A simulation program has been written here which produces 
the correct waveform incorporating the sampled point- target response effects and 
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fits the "Incorrect but fast" model waveform to these points. By varying the 
Input parameter values systematically and comparing the results, the desired 
correction curves can be obtained. Figures 9 through 12 show the needed correc- 
tion curves for SWH, time origin ("track point"), off-nadir attitude, and 
surface skewness, respectively, for 63-po1nt, six-parameter, uniformly-weighted 
data fitting with the skewness and attitude angle lightly constrained. Any 
change In these conditions, number of points or parameters fitted and so forth, 
would necessitate rerunning the entire simulation to generate a new set of 
correction curves for the new conditions. There 1r> nothing In the generation of 
the correction curves to account for errors In Individual waveform sampler 
gains; that possible problem must already have been solved. A built-in 
assumption In generating Figures 9 through 12 was that for a given set of Input 
parameters, the average of a number of waveform fits from waveform samples with 
typical random noise would converge to the zero-noise result. Figures 9 through 
12 were generated for zero-noise waveform samples, equivalent to Infinite 
averaging time, but for the ten second, 1000 pulse averages of actual SEASAT-1 
data that we have been using the noise standard deviation Is around 3% of the 
mean level In each sampler. We have since added realistic random noise to the 
simulation program and the mean values of the fitted parameters from the noisy 
simulated data do seem to converge to the results of the zero-noise case but 
this conclusion is based on a limited number of runs and more work Is needed on 
this question. 


Figure 9. SEASAT-1 SWH corrections. 
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WAVEFORM SAMPLER GAIN ADJUSTMENTS 


Preliminary results from our first waveform fitting to the SEASAT-1 data 
Indicated that there were Individual sampling data biases and that some method 
had to be developed for "fine-tuning" the Individual sampler gains after using 
the Calibration Mode II data as a starting point for these gains. Obvious 
patterns were visible In the residuals after the waveform fitting; some of the 
sampler's residuals appeared always to be positive, others negative, and we 
Intended using averages of these residuals over a number of sets of Input data 
to fine-tune the starting gains. This situation Is similar to the earlier 
GEOS-3 radar altimeter, although there additive gain biases were assumed rather 
than the present cases' multiplicative gain biases because of the hardware 
differences. Walsh, for example, described a GEOS-3 sampler bias adjustment 
procedure in Ref. 5, 

Figure 13 is a summary of the gain bias fine-tuning which we Initially 
Intended using. This procedure was In fact used In our data reduction just prior 
to the GOASEX II Workshop In June 1979. and In the Wallops contribution to the 
altimeter portion of the report from the workshop (Ref. 3). An implicit assump- 
tion for the process In Figure 13 Is that the gain bias adjustments needed are 
relatively small and that the Individual sampler final gains are distributed 
more or less randomly around the starting-point gains from Calibration Mode II. 
There is also the requirement that the off-nadir angle be known, and for the 
early June 1979 work we used the SEASAT attitude control system-derived angle as 
available from the JPL-processed data tapes. Figure 13 Indicates an optional 
final adjustment over the full sampler set, but we have not used that option In 
work to date. 
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Figure 13. Summary of gain fine-tuning procedure. 
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I. GAIN CORRECTIONS OUTSIDE CENTRAL PORTION OF SAMPLER SET 

a. Select Low SWH cases only from data set 

b. For known attitude, fit waveform function 

c. Use residuals outside region A-B as multiplicative gain corrections 

II. GAIN CORRECTIONS IN CENTRAL PORTION OF SAMPLER SET 

a. Select High SWH cases only from data set 

b. For known attitude, fit waveform function 

c. Use residuals inside region A-B as multiplicative gain corrections 

HI. (OPTIONAL) FINAL ADJUSTMENT OVER ENTIRE SAMPLER SET 

a. For known attitude, fit waveform function 

b. Use residuals over all samplers as final multiplicative corrections 
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Figure 14a shows the starting-point sampler gains from an average of Cali- 
bration Mode II data, and Figure 14b shows the final set of gains after a pro- 
cedure similar to that outlined In Figure 13. Notice the large (and, to this 
date, unexplained) changes In the gains particularly In the later plateau por- 
tion. For discussion later In the section on a plateau-region attitude estima- 
tor, notice the (negative) non-zero slope In the Figure 14a smoothed sampler 
gain vs sampler number curve In the region of waveform samplers 39-56. The 
corresponding slope In Figure 14b Is zero or slightly positive. Notice also 
that the final gain set tends to be smoother gate-to-gate than the starting- 
point set. The two sets of gains In Figure 14 are so different that we might 
almost as well have started by assigning unit gains to all samplers as the 
starting point for carrying out the gain fine-tuning. In other words, the 
Calibration Mode II data are almost useless to us. 

AHITUDE ANGLE EFFECTS 

When we applied the final sampler gain set described above to a variety of 
experimental data in the work In early summer 1979, we found a number of regions 
of over-ocean SEASAT-1 data in which the waveform fit results were poor because 

the sampled plateau region decayed faster than allowed even at 0° off-nadir 

> 

angle (recall that the fastest decay is at 0® as shown in Figure 4). The six- 
parameter waveform fitting process could not reach reliable estimates of the 
other parameters because of the difficulties in the fitted off-nadir angle in 
this situation. We were also finding that the off-nadir angle versus time 
results, even in regions where the fit was not trying to produce negative 
absolute values, did not correlate well with the attitude versus time results 
from the attitude control system's attitude estimates. These various observa- 
tions cast considerable doubt on our ability to use the attitude control system's 
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Figure 14a. Initial Sampler Gains from Calibration Mode II 

(Solid line is 5-point smoother with weights 1:3: 5:3:1) 
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SEASAT Waveform Sampler Number 


attitude estimates for input to the gain fine-tuning procedure. Figures 17-19, 
in the following report section describing a possible plateau-region attitude 
estimator, show our recent results supporting our claims of attitude estimation 
difficulties. We emphasize here that we need reliable attitude estimates only 
for the initial process of fine-tuning the individual waveform sampler gains. 
Once these on-orbit opera te-status gains are determined, we assume that they are 
time-stable and we no longer require attitude estimates because the attitude can 
be one of the parameters fitted in the rest of the waveform processing. 

Figure 15 is a schematic diagram to accompany this paragraph's discussion 
of interaction of the gain uncertainties and the attitude uncertainties. 

Suppose for example that only sample points at times earlier than B (Region I 
only in Figure 15) are fitted by a several -parameter waveform and that the 
SEASAT-1 attitude control system-estimated attitude’ 5^ is used as a fixed, 
known attitude value. Some sort of gain adjustment over all waveform samples 
will then, in effect, force the plateau region (Regions II and III in Figure) 
waveform samples to lie on the curve labeled in this figure. Suppose however 
that the true attitude value were instead, where > Sg* Ih this case the 
plateau-region final gains would all be in error in such a way that when the 
attitude angle is treated as a fitting parameter the fitted attitude values 
would all be lower than the true attitude values. In data spans where the true 
attitude itself is near zero, the waveform fitting procedure would then try to 
attain the impossible case of a negative absolute attitude value, and all other 
fit parameters would then be in error also. The situation just described is 
what we think did happen in the work in the early summer of 1979. 
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Waveform Amplitude 


Figure 15. Figure for Discussion of Gain Biases and 
Pointing Angle Relationship. 



PLATEAU-REGION ATTITUDE ESTIMATOR 


We are in the difficult position of believing that a sampler gain adjust- 
ment procedure is required and yet not having a reliable enough attitude knowl- 
edge to carry out the gain adjustment procedure. Once the gains are adjusted, 
our own waveform fitting can then vary the attitude as one of the fitted param- 
eters but not before this adjustment has been completed. 

The SEASAT-1 altimeter had a wide sampling gate later on the waveform 
plateau. This attitude gate, in a waveform region sensitive to attitude angle, 
was intended to provide data for attitude estimation. There was a similar 
attitude gate on GEOS-3 (sometimes referred to as the "attitude/specular gate" 
in the GEOS-3 literature) and these gate's locations are indicated on Figure 1. 
Unfortunately the SEASAT-1 attitude gate's scale factors were incorrectly set 
and the resulting attitude gate data are always at or near saturation. The 
result is that the attitude gates output is insensitive to attitude changes and 
thus is not useful to us as an attitude estimator. 

The following several paragraphs sketch a "quick and dirty" attitude esti- 
mator using plateau region waveform sampler data. The estimator has calibration 
difficulties of its own but it should provide us with at least a look at the 
attitude trend and to identify data expanses in which the attitude is relatively 
constant. These constant attitude regions will then be input data for a revised 
sampler gain setting procedure. The basic procedure is to fit an exp(-6t) 
functional form to plateau region waveform samplers, and to estimate the atti- 
tude from simulation generated curves of 6 vs 

Recall that, as shown in Figure 1, the flat-sea impulse response function 
for the nadir-pointing (C = 0) SEASAT-1 radar altimeter has an exponential decay 
form at times greater than zero. For non-zero ?, the behavior is not a simple 
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exponential but can still be fitted in a local region by an exponential. We 
allow for a noise baseline and assume that for low to moderate SWH the Region II 
waveform in Figure 15 may be represented by 

W(t) * A exp(-6t) + b , 

We form the variable Z by 

Z = tn W(t) - b , 

and we do a least-squares fit of data from waveform samplers 39 through 56 to 
determine 5 from 

Z = tn A - «t . 

The b is estimated from the average of samplers 4 through 13, and we usually use 

a data smoothing interval of 10 seconds for each 6 estimated. 

For low to moderate SWH, the region of samplers 39-56 is well away from 

where the point-target response shape details affect the shape of the mean 

waveform, so there is no real difficulty caused by using the Ref. 9 expressions 

as a complete waveform description for this plateau region. In a separate 

simulation program we have generated an expected return waveform and then fitted 

the above exponential function to generate the curves of 6 vs ^ given in Figure 16. 

The 6 estimated from experimental data will depend on what sampler gains 

are used for the samplers 39-56. We used the five-point-smoothed gains from 

Figure 14a. From the calculated Pp 5 (t) at nadir point ($ = 0) as in Figure 1, 

-3 -1 

we know that the theoretical maximum value for 6 is = 2.66x10 ns . 

m 
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Attltude-Estlmatlrig FufiCtlon, 10^x6, in ns 


Figure 16. Attitude-estimating function vs attitude 
angle as function of SWH. 
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There were many regions of over-ocean SEASAT-1 data for which the experimental 6 

values were greater than and we believe this was due to the operate-only 

gain trend (that is, to the different trend in the sampler gain vs sampler 

number in samplers 39-56 for Figure 14b as compared to Figure 14a). The first- 

order effect of a built-in erroneous gain vs sampler number trend would be a 

shift in the resulting fitted 5 values. The various data we have looked through 

-3 

indicate that the maximum experimental 6 values are probably around 3.6x10 ns 

-3 -1 

so we subtract 0.94x10 ns from all experimental £ before entering the curves 
of Figure 16 with 6 and SWH to get an estimate of 5. The first-order result of 
this 0.94 subtraction is a shift in the estimated C values. 

Figures 17-19 show smoothed 5-est1mated attitudes together with the atti- 
tude derived from the attitude control system. Figure 17 shows the two atti- 
tudes in fairly reasonable agreement in the first half of the interval shown, 
but with local maximum in different locations in the second part of the plotted 
interval. Figure 18 shows a different pass with more disagreement in locations 
of local maxima and minima. In neither Figure 17 nor Figure 18 is there strong 
indication of a bias or offset between the two different attitude estimates; 
such an offset would be the first effect of an error in the subtractive 6 
correction. Figure 19 however does show such a relative shift or offset; and 
the differences shown on Figure 19, particularly around 345 minutes, are signifi 
cant to our waveform processing. 

In summary, the saturation of the SEASAT-1 attitude gate prevents its use 
for attitude estimation. We have described the way in which an attitude esti- 
mator can be developed using plateau region sampler values (plus low-numbered 
samplers to estimate the baseline). This attitude estimator is not absolutely 
calibrated because of sampler gain problems, but should nonetheless be useful 
for attitude trend measurements. 


time 



Time, minutes 




Figure 19. SEASAT-1 estimates of off-nadir angle vs time. 
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SUMMARY 


Our work has been directed toward fitting all 63 SEASAT sampling gate 
values to a six-parameter model waveform function, once a process of adjusting 
the Individual sampler gains has been carried out. Ideally, this gain "fine- 
tuning” would not be necessary; practically, we think it is very Important for 
the SEASAT situation. The model waveform for a nearly-Gausslan system point 
target response has been described elsewhere (Ref. 9) together with a fuller 
discussion of the background and assumptions for the radar altimeter mean return 
waveform calculations. As shown by Figure 6 the SEASAT system point target 
response Is not describable as nearly-Gausslan, and the actual sampled polnt- 
taro'it response must be convolved with the general functional form from Ref. 9 
(representing "everything else," everything except the point-target response) to 
obtain the complete model waveform to be fitted to the 63 SEASAT waveform 
samples. The waveform fitting procedure itself is completely straightforward; 

It was sketched In Ref. 7 for a four-parameter situation and the extension to 
six parameters is trivial. Our Fortran source program for the waveform fitting 
contains a variety of minor adjustments which specifically adapt a once-general 
fitting program to the special requirements of the SEASAT case, but these 
changes are generally indicated by comments in our source listing (available now 
on request, and to be published soon as a separate report). 

The considerable difficulties which we have encountered have been in three 
general areas; 1) obtaining at final individual waveform sampler gains and 
understanding why those gains differ from initial values from Calibration 
Mode II, 11) incorporating the actual system point- target response into the 
waveform model (and the knowledge of the point- target response also depends on 
resolving sampler gain questions) and iii) obtaining reliable values for the 
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off-nad1r attitude angle, at least for the initial sampler gain-adjusting phase 
of the work. Once sampler questions are resolved, incorporating the detailed 
system point-target response function into the waveform description is simple 
enough with the main difficulty being the practical one of markedly increased 
computation time resulting from a numerical convolution step within the fitting 
procedure. The question of individual waveform sampler biases remaining after 
using the best available calibration data is not unique to SEASAT; we had a 
similar situation in previous GEOS-3 work (where the waveform sampler biases 
were taken as additive rather than multiplicative) and, earlier, in work with 
the Skylab S-193 radar altimeter. Clearly if there were no difficulties indi- 
cated with the sampler individual gains and consequently no "fine-tuning" pro- 
cedure necessary, then we would not require that the off-nadir attitude angle be 
known because this could be found as one of the fitted waveform parameters. 

We would hope that our SEASAT-1 work to date would underscore the impor- 
tance of good design of a relatively clean transmitted pulse with very dense 
sampling of the shape of the system point-target response, and would emphasize 
the need for great care in the design of the on-orbit calibration steps. The 
loss of useful SEASAT-1 attitude gate data points to one specific area for 
improvement in future systems. We have had difficulty using the on-orbit 
calibration data to interpret altimeter's waveform data while tracking (notably 
in the relative waveform sampler gains of Calibration Mode II vs the apparent 
gain fall-off in the late plateau samplers). Additional flexibility in the 
tracking operation, such as the ability to change the tracking point within the 
waveform sampler set, could be very useful for future radar altimeters. 
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